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ABSTRACT 

We present three Suzaku observations of the black hole candidate 4U 1957+11 (V1408 Aql) — a source 
that exhibits some of the simplest and cleanest examples of soft, disk-dominated spectra. 4U 1957+11 also 
presents among the highest peak temperatures found from disk-dominated spectra. Such temperatures may be 
associated with rapid black hole spin. The 4U 1957+11 spectra also require a very low normalization, which 
can be explained by a combination of small inner disk radius and a large distance (> lOkpc) which places 
4U 1957+11 well into the Galactic halo. We perform joint fits to the Suzaku spectra with both relativistic 
and Comptonized disk models. Assuming a low mass black hole and the nearest distance (3Mq, lOkpc), the 
dimensionless spin parameter a* = Jc/GM^ > 0.9. Higher masses and farther distances yield a* « 1. Similar 
conclusions are reached with Comptonization models; they imply a combination of small inner disk radii (or, 
equivalently, rapid spin) and large distance. Low spin cannot be recovered unless 4U 1957+11 is a low mass 
black hole that is at the unusually large distance of > 40kpc. We speculate whether the suggested maximal 
spin is related to how the system came to reside in the halo. 

Subject headings: accretion, accretion disks - black hole physics - radiation mechanisms:thermal - X- 
rays;binaries 



1. INTRODUCTION 

4U 1957+11 is one of the few black hole candidates 
(BHC), like LMC X-1, LMC X-3, and Cyg X-1, that histor- 
ically has been persistently active. In contrast to these other 
sources, however, it i s a Low Mass X-ray Binary (LMXB ; 
see Thorstensen 198i lHakala et al.l[T999l:lRussell et aUlMol: 
[Bayless et al. 2011). Also in contrast to other persistent 
Galactic BHC sources, 4U 1957 +11 apparently h as remained 
always in a spectrally soft state (N owak et al . 2008). Its soft 
spectrum is well-modeled as a simple disk (H 
IWijnands et all llool iNowak et alJ 12008 
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a multi-temperature blackbody (.Mitsuda et al 



1984) characterized by a peak temperature and a normaliza- 
tion related to the disk inner radius, inclination, and object 
distance. Approximately 15% of the over 40 pointed obser- 
vations taken with the i?ossi X-ray Timing Explorer (RXTE), 
however, also show a significant contribution (20-75% of the 
observed 3- 18keV flux) from a non-thermal/stee p powerlaw 
component (No wak et alj|200l iDuimet al.ll2010t) . 

Optical observations suggest that we may be viewing the 
disk in 4U 1957+11 at a high incHnation of - 75°. Mod- 
ulation over a 9 .33 hr orbital period has ranged from ±10% 
and sinusoidal (lTl i orstense"nl il987) to ±30% and complex 
(lHakalaetal.1 119991) . The latter has been interpreted as a 
high inclination warped disk b eing partly oc c ulted by the sec- 
ondary. On the other hand, iRussell et al.l (12010 ) failed to 
find any opti c al mo dulation over the orbital period, while 
iBavless et al.l ( 1201 ll) found nearly sinusoidal modulation. 
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These latter authors attributed the optical modulation to illu- 
mination of the secondary, and noting a lack of any flattening 
at the lightcurve minimum, they modeled the system inclina- 
tion with values as low as 20° and not greater than 70°. 

Currently, the mass, distance, and inclination of 
4U 1957+11 are unknown. Based upon a comparison 
of the optical flux to the estimated optic al luminosities of 
other BHC, however, IRussell et al.l ( 120101) have argued that 
4U 1957+1 1 must lie at a distance > 20kpc. The possibility 
of a high in clination, however, makes such comparisons 
problematic. IRussell et alj (|201 V) recently have presented 
near-simultaneous radio (Expanded Very Large Array - 
EVLA), optical (Faulkes Telescope North), and UV/X-ray 
(Swift) observations of 4U 1957+11. These observations 
show the soft state flux to be extremely quenched (factors of 
330-810 times lower than the expected radio/X-ray flux ratio 
for hard state black hole sources), and demonstrate that the 
optical/UV flux are roughly consistent with the expectation 
of an irradiated disk. Those observations, however, do not 
definitively constrain the source's inclination or distance. A 
minimum distance of > 5 kpc; however, has been determined 
via high resolution X-ray spectroscopy. Ne IX 13 .45 A 
absorption, associated with the warm/hot phase of the 
interstellar medium, is detected with a sufficient equivalent 
width to suggest that 4U 1957+11 resides a bove the galactic 
plane at this minimum distance or beyond dYao et al.l 120081 : 
Nowak et al. 2008). 

INowak et alj ((2008) showed that, with the exception of the 
« 15% of the time when a non-thermal component becomes 
prominent, disk models of the soft X-ray component have 
an approximately constant radiu s but a temperature that in- 
creases with flux (see Fig. 9 of INowak et alj 120081) . There 
is some indication for a slight increase in disk radius at the 
lowest flux levels, which was interpreted as possibly being 
associated with the source approaching a transition to a spec- 
trally hard state. Such a transition is expected to occur at 
2-3%LEdd, where LEdd is the Eddington luminosity for the 
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source in question ( lMaccaroli3l2003[) . At the other extremes 
of flux in 4U 1957+1 1, the simple disk spectrum "drops out" 
(disk radius decreases in disk+powerlaw fits, or disk temper- 
ature decreases in disk+Compton corona fits) and is replaced 
by the non-thermal/steep powerlaw component dNowak et alJ 
12)08). 

Although other persistent systems transiently show disk- 
dominated spectr^ the persistent, soft spectra of 4U 1957+1 1 
are perhaps the simplest and cleanest examples of "disk spec- 
tra". XMM-Newton and Chandra observations show that 
there is very little absorption (A^h = 1-2 x 10^' cm"^), while 
RXTE observations show that a hard tail contributes usually 
< 20% of the flux. Furthermore, any additional components 
such as a broad Fe line (typically < 2% peak residual in the 
fits) or a "smeared Fe edge" (which was never required o r in- 
dicated in any of the RXTE spectra; iNowak et alj|2008b are 
very weak or absent. Thus the soft spectrum of 4U 1957+1 1 
becomes the ideal testbed for modern disk atmosphere mod- 
els that incorporate spin and other General Relativistic effects 
into their calculations_(e.g.. lLi et"ani2005l; iDavis et al.ll2005L 
|2006; Shafee et al. 2(mt). 

[Nowak et al.i (i2008ir showed that the disk-fits to the XMM- 
Newton, Chandra, and RXTE spectra of 4U 1957+11 are 
characterized by a very low normalization, indicating some 
combination of large distance and low compact object mass, 
and very high inner disk temperature (1.3-1. 8 keV). In other 
sources the latter has been associated with high black hole 
spin. To assess this possibility, Nowak etal. (2008) applied 
the relativistic disk model, kerrbb, of Li et al. (2005). The 
model parameters include a compact object mass and dis- 
tance, disk inclination, accretion rate, spectral hardening fac- 
tor (the ratio of color temperature to effective temperature), 
and dimensionless spin of the black hole, a* = Jc/GM-. For 
a wide variety of masses and distances, fits to the XMM- 
Newton and Chandra spectra preferred maximal spin, a* 
The RXTE spectra, in fact, could not be described with low 
spin, even when including a coronal or powerlaw component. 

Various theoretical models of a disk atmosphere have sug- 
gested a spectral hardening factor of ^ 1.7. If this is correct, 
and utilizing the fact that the lowest flux observations must 
be > 3%LEdd, then the spectral fits of 4U 1957+11 strongly 
suggest that it is a maxi mally spinning, sal 6M(7i black hole at 
a distance of «22kpc (INowak et al.ll2008h . In this work we 
explore that claim with three recent Suzaku observations. 

In Section |2] we describe the Suzaku observations and the 
details of the data processing. We then present the model 
fits to the resulting spectra in Section [3] Section |4] describes 
the scaling relations that we used to translate fits using a low 
mass, short distance black hole to comparably good fits in- 
volving a larger mass, further distance black hole. We then 
present our conclusions in Section|5] 

2. OBSERVATIONS AND DATA ANALYSIS 

The Suzaku data were reduced with tools from the HEA- 
SOFT v6.9 package and the calibration fil es current as of 20 1 
September. The instruments on Suzaku (Mitsuda et al. 2007) 
are th e X-ray Imaging Spectrometer (XIS; Koyamaetal. 
|2007^ CCD detector coverin g the » 0.3-lOkeV ba nd, and the 
Hard X-ray Detector (HXD: lTkkahashi et al.l200'7h comprised 
of the PIN diode detector (PIN) covering the « 10-70keV 

' LMC X- 1 is likely wind-f ed, and, for a soft s tate, shows an unusual and 
highly variable disk spectram (Nowak et al.'2001'). LMC X-3 cycles thi'ough 
soft/hard state transitions (Wilms et al. 2001). 



TABLE 1 
Log of 4U 1957+1 1 Observations 



Date 




ObsID 


Exposure 


(yyyy-mm-dd) 


(MJD) 




(ksec) 


2010-05-04 


55320.82 


405057010 


35.80 


2010-05-17 


55333.80 


405057020 


37.20 


2010-11-01 


55502.19 


405057030 


15.50 



Note. — Exposure times are after good time filter- 
ing, and represent the maximum exposure in an XIS 
detector when combining all data modes. MJD is 
the (non-barycenter corrected) average Modified Juhan 
Date of the observation. 




Fig. 1. — The aspect corrected XIS 1 0.5-8 kev image on a logarithmic 
intensity scale from ObsID 405057010. The boxes outline the detector re- 
gion from which the spectra were generated. The length of the short edge of 
the outer box is approximately 4'. 5. The inner boxed regions with diagonal 
slashes were excluded from the spectra. 

band and the gadolinium silicate crystal detector (GSO) cov- 
ering the « 60-600 keV band. The XIS had four separate 
detectors, XIS 0-3, with XIS 1 being a backside illuminated 
CCD. XIS 2 was lost due to a micrometeor hit in late 2006, 
and thus we only consider data from the remaining three XIS 
detectors. For all three observations described here, the hard 
X-ray fluxes were too low to lead to a significant detection in 
either of the HXD detectors, and therefore we do not include 
these data. A log of the observations is found in Table [T] 

All of the XIS observations were run in a sub-array mode 
wherein 1/4 of the CCD was exposed with 2 second inte- 
gration times. In preparing the XIS spectra, we first cor- 
rected each detector for Charge Transfer Inefficiency (CTI) 
using the xispi tool, and then reprocessed the data with 
xselect using the standard xisrepro selection criteria. 
Due to thermal flexing of the spacecraft, the attitude of the 
Suzaku spacecraft exhibits variability over the course of the 
observations and therefore the image of the source is not at a 
fixed position on the CCD. Standard processing reduces this 
variab ility and improves th e Point Spread Function (PSF) im- 
age (lUchivama et al.l 12008 ): however, the standard tools do 
not yet fully correct the variable attitude induced blurring. We 
created an improved attitude reconstruction, and hence nar- 
rower PSF images, using the aeattcor . si software de- 



Suzaku Observations of 4U 1957+1 1 



3 




0.5 1 2 5 

Energy (keV) 



Fig. 2. — The combined XIS 0, 1, and 3 spectra from each of the three 
Suzaku observations of 4U 1957+11, flux corrected with the detector re- 
sponses. (ObsID 405057010 is the top spectrum, obsID 405057020 is the 
middle spectrum, and obsID 405057030 is the bottom spectrum.) The three 
observations share a common binning of approximately half width half max- 
imum resolution of the detectors. Lighter colored data points represent 
energy bands that were ignored in the spectral fits described below. 

scribec0 by Nowak et al.l ( 1201 1). For the third observation, 
however, approximately half of the good time intervals (i.e., 
half of each spacecraft orbit) were affected by wobbles severe 
enough (> 1') such that a large fraction of the image on the 
XIS 1 detector fell beyond the boundaries of the detector's 
active regions. Although the images did not fall beyond the 
active region boundaries on the XIS and XIS 3 detectors, 
there are serious questions about the accuracy of the gener- 
ated response functions for such large off axis angles. We 
therefore excluded these times from the spectra, and reflect 
this in Table [U 

Using the image generated from the improved attitude cor- 
rection, w£_flien_usedjhe pile_estimate . si S-Lang 
scripl3 dNowak et alj|201 fl) to estimate the degree of pileup in 
the XIS observations. For the spectra described in this work, 
the center of the PSF images were affected by as much as a 
33% pileup fraction. We therefore exclude two, overlapping 
rectangular regions (« 15 x 60 pixels, nearly at right angles to 
one another) aligned along the brightest portions of the cen- 
tral PSF (see Fig. [Til. We estimate that the remaining residual 
pileup fractions are < 3%. 

The outer limit of the extraction area was a rectangular re- 
gion limited by the 4' width of the 1/4 array boundaries. Ad- 
ditionally, we extracted background spectra from small rect- 
angular regions, near the edges of the chips, that avoided as 
much as possible both the 4U 1957+11 image as well as the 
detector images of the on-board calibration sources. It is not 
entirely clear to what degree these regions are dominated by 
"true background", as opposed to the wings of the PSF from 
the source. However, the normalized background flux remains 
at < 2% of the source flux between 0.5-8 keV, and thus does 
not strongly affect the analyses presented in this work. The 
background flux starts to become an increasingly larger frac- 
tion of the source flux above 8 keV. For this reason, as well as 
due to any other uncertainties in the high energy calibration 
of the XIS detectors, we do not consider data above 8 keV. 

Events in the XIS detectors were read out from either 3x3 
or 5 X 5 pixel islands. We created individual spectra and re- 
sponse files for each detector and data mode combination. Lit- 

^ See also: jhttp : / /space .mit ■ edu/CXC/ software/ suzaku/| 
aeatt . html. 

^ |http : //space .mit ■ edu/CXC/software/suzaku/pest ■ htinll 




Energy (keV) 

Fig. 3. — Ratios of the flux corrected spectra shown in Fig.|2] The top ratio 
is observation 1/observation 3, the middle line is observation l/observation2, 
and the bottom line is observation 2/observation 3. 

tie source variability was detected over the course of each ob- 
servation, therefore the spectra were created for the full good 
time intervals represented by standard processing (excluding 
the times of very bad attitude errors, discussed above). Re- 
sponse matrices and effective area files were created with the 
xisrmfgen and xissimarfgen tools, respectively. For 
each observation, spectra for the different readout modes were 
combined during analysis on a chip-by-chip basis using the 
combine_datasets function within the Interactive 
Spectr al Interpretation System (ISIS) fitting 
program (Houck & Denicola' '2000). That is, observations 
were kept separate, but within an observation all XIS spec- 
tra were combined together, all XIS 1 spectra were combined, 
etc. 

Although spectra for each observation and each XIS chip 
were kept separate, all spectra were jointly grouped on a 
common grid such that the XIS spectra for the third (i.e., 
the faintest and shortest integration time) observation had a 
minimum combined signal-to-noise ratio of 5 in each energy 
bin (p» 25 total counts in each bin, although the estimated 
background was included in the S/N criterion) and that the 
minimum number of channels per energy bin was at least 
the half width half maximum of the spectral resolution (see 
iNowak et al. 201 1. for a further explanation of this criterion). 

In Fig. |2] we present the flux corrected spectra, i.e., the 
counts spectra divided by the integrated response (see the 
ISIS user manual for a full description of this process). The 
spectra show strong deviations from a smooth continuum in 
the 1.5-2.6keV region, near X-ray spectral features of Si and 
Ir. These systematic effects are seen in a wide variety of un- 
related Suzaku spectra, hence our decision to exclude them 
from our fits. Likewise, there is a strong spectral feature near 
0.55 keV, which we also suspect of being systematic in na- 
ture, hence our exclusion of data below 0.7 keV. To avoid 
these regions of poorly understood response, we only consid- 
ered spectral energy ranges 0.7-1.5 keV and 2.6-8 keV. The 
smaller features near 0.75 keV and 0.9 keV are plausibly as- 
sociated with Fe L, Ne K, and Ne IX absorption in the inter- 
stellar medium (ISM), as we discuss below. 

The suspected systematic features are relatively constant 
among the three observations, as is evident from the ratio of 
the flux corrected spectra as shown in Fig. [3] None of the 
systematic features are seen in these ratios. The ratios are rel- 
atively flat < 2 keV, and show a curve upward, with the degree 
of this curving being greatest for the ratio of the flux corrected 
spectra for the first and third observation. We discuss this fur- 
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Fig. 4. — Expectations for the flux ratio in a diskpn model when solely 
changing the peak temperature of the disk. The reference model has a 1.3 keV 
peak temperature. The peak temperature is varied in steps of 0.05 keV, up to 
1.5 keV. 




Energy (keV) 

Fig. 5. — The joint fit of an eqpair Comptonization model to the three 
spectra. The spectra share a common interstellar absorption column and a 
common normalization (i.e., disk inner radius), but have independent disk 
temperatures and coronal pai'ameters. Black/grey spectra are the combined 
XIS and 3 spectra, and blue/light blue are the XIS 1 spectra. Lightly shaded 
data were not included in the fit. The top residuals panel shows the ratio resid- 
uals for obsID 405057010, the middle residuals panel shows the ratio residu- 
als for obsID 405057020, and the bottom panel shows the ratio residuals for 
obsID 405057030. 

ther in the next section. 

3. SPECTRAL FITS 

The spectral shape shown in Fig. |2] and the ratios seen 
in Fig. [3] are those expected for a disk-dominated spectrum 
wherein the flux changes are primarily dominated by vari- 
ations of the peak temperature of the disk. We show the 
theoretical expectation fo r such a disk in Fig. |4] Here we 
use the diskpn model (iGierliiiski et alJ IT999) . which es- 
sentially is the spectrum f rom a Shakura-Sunyaev type disk 
(IShakura & SunyaevI 1973 1) calculated using a simple pseudo- 
Netwonian potential and a no-torque boundary condition at 
the disk inner edge located at 6GM/c^. The flux goes to 
zero at the disk inner edge, and the disk temperature peaks at 
w IQGMc^. This basic disk spectrum will underly the Comp- 
tonization fits discussed below. 

At energies approximately below the peak temperature, 




4 5 6 7 8 

Energy (keV) 



Fig. 6. — A simultaneous fit to all three observations using the eqpair 
model, but excluding the broad Gaussian line in the Fe region. (The spectra 
were combined for the plot, but not the fits. 

dominated by the Rayleigh- Jeans part of the thermal spec- 
tra, the ratios are fairly flat. At energies above the disk peak 
temperature, on the Wien tail of the thermal spectra, the ratios 
curve upward. There is a very good correspondence between 
the curves in Figs. [3] and |4] However, based upon our ex- 
perience with RXTE spectra which extended out to « 20 ke V, 
there is the possibility of additional spectral hardening beyond 
that associated with any rise in disk temperature. We therefore 
consider the eqpair Comptonization model (Coppi i 19991) 
which can use the di skpn as its input seed photon spectrum. 

Given that the three spectra appear commensurate with a 
constant disk radius, and that such a nearly constant disk 
radius is consistent with our experience for w 85% of the 
RXTE spectra that we have previously examined (see Fig. 9 
of Nowak etal. 2008), we explore a model wherein we fit 
all three spectra simultaneously with an eqpair model for 
which we constrain the eqpair normalization (i.e., the inner 
disk radius) to be the same for each spectrum. The physical 
implications of this common normalization are discussed in 

m 

In addition to constraining the eqpair normalization, we 
also constrain the absorption to be the same for each obser- 
vation, which we here model w ith the TBnew mo del (an up- 
dated version of the model of I Wilms et alJ l2000l using the 
interstellar abundances from that work) and using a fixed Ne 
IX 13.447 A absorption line with equivalent width of 7 A. The 
latter is co nsistent wi th the high sp ectral resolution studies of 
lYao et alj (|2008) and iNowak et alJ (2008) ( see the below dis- 
cussion of the ISM absorption). Excluding this component 
increases the fitted by « 12. 

For each observation we allow the disk peak temperature, 
fcrdisk, the coronal compactness, £h/£s, and the coronal seed 
optical deptlQ, r,, to be independent parameters. We include 
Gaussian lines to the fits; however, as their presence is not 
strongly confirmed and there is no good evidence for a nar- 
row line in either these spectra nor in the previously exam- 
ined Chandra -HETG spectra ( Nowak et al...2008) , we tie all 
the line energies (constrained to lie between 6.2-6.9 keV) and 
widths (fixed to 0.3 keV), and let only the line normalizations 

^ Table[2] however, presents the net optical depth, as eqpair includes 
calculation of the equilibrium pair production in the corona. 
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TABLE 2 

constant*TBnew* (eqpair + gaussian) FITS TO 4U 1957+1 1 SPECTRA 



ObsID 


(lO^l 


Weqp 
(10^) 


(keV) 


4/4 


Tnct 


N-, 
(10^) 


QisO 


Cxi.s3 


0.7-8 keV Flux 
(10"' ergcm"^ s) 


xVDoF 


405057010 
405057020 
405057030 


200^ °°' 
"■^""-0.001 


^•^^"-0.001 


, 477+O.001 
'-0.001 
1 ^42+0.001 

1 310+°"'" 


37+<'-''^^ 
"•^'-0.02 

05+O '" 
"•"-^-0.01 

04+° '" 
"•"^-0.01 


Q Qiy+fl.002 
^'■^'^ '-0.002 
2 4+0.1 

44+0.1 


1 4+0.4 

l+°-2 

1 3+°-'* 
^•^-0.4 


1 019+0.002 
^■"^ -0.002 
1 078+0.002 
^■'''°-0.002 

'■0°'-0.002 


1 020+°°°2 
^•"^ -0.002 

1 nR^+0.002 
1 057+*°°' 


1.09±0.01 
0.80 ±0.04 
0.72 ±0.03 


1750.4/1014 


405057010 
405057020 
405057030 


-0.008 


1 926+0.002 

^•^^"-0.002 


, 476+0.001 
^■^'"-0.001 

I ,41+0.001 
1 310+°"'" 


A ,0+0.03 

0.3S_(,o3 

05+° 

05+° 
"•"-^-0.01 


Q Ajy+0.002 
'-0.002 
2 4+0.1 

45+0.1 


1 4+0.4 

1+°-'' 
"■'-0.1 

, ,+0.4 
^•-'-0.4 


1 020+°''°2 
1 077+0002 
-0.002 

1 066+° °°^ 

^■"°°-0.004 


1 020+<'°''3 

1 riRl+0 003 

^•"'*'-o.ooi 
1 056+*°°'* 




1016.9/1014 



Note. — EiTors ai'e 90% confidence level for one interesting parameter, except for the 0.8-8 keV (absorbed) flux error which is from the standard deviation 
of the XIS detector cross normalizations. Nh is the neutral column; A'cqp is the eqpair normaHzation (see text); kTusi^ is the peak temperature of the disk seed 
photons; £h/^s is the relative corona to disk compactness (£s is fixed to 1); Tnei is the net coronal optical depth (including pair production); Nry is the Gaussian 
line normalization in units of photons cm"^ s"' ; and CxisN are the normalization constants for the XIS and 3 detectors, relative to the XIS I detector The 
first group of parameters are for data that include only statistical errors, while the second group of parameters are for data that also include 1.35% systematic 
enors. 



remain free. 

The results of these fits are presented in Fig. |5] and the pa- 
rameters are given in Table|2] Overall, the quality of the fits is 
very good, yielding xl/i' = 1756/1014= 1.74. Although the 
reduced is not 1, the addition of only 1.35% systematic un- 
certainties to each dataset brings the value down to that level, 
without significantly altering the best fit parameters, as also 
shown in Table |2] The addition of systematic uncertainties 
slightly increases the parameter error bars. 

Given the obvious systematic differences between the 
frontside (XIS and 3) and backside (XIS 1) detectors, seen 
in the residuals presented in Fig.|5] we speculate that our fits 
are near "optimal' . Throughout the rest of this work, the fig- 
ures will include only statistical errors; however, tables will 
include sets of parameters with error bars that were calcu- 
lated for spectra to which 1.35% systematic errors have been 
applied. This systematic error amplitude was chosen on the 
basis that it yields a reduced value of w 1 for the fits de- 
scribed above, and we retain this level in subsequent model 
fits for ease of comparison. 

The primary difference among the three spectra are the disk 
peak temperatures, which range between 1.31-1.48keV. The 
coronal contribution is relatively weak, given the low values 
of the coronal compactness compared to the disk compact- 
ness, ih/is- It is interesting to note that the two fainter ob- 
servations indicate less energetic coronae — ^h/^s values are 
lower — but higher optical depth ones. (Coronal compact- 
ness is its input energy divided by a characteristic radius, here 
measured relative to a compactness for the disk. If geomet- 
ric changes are absent, then a lower fitted compactness is in- 
dicating a less energetic corona.) The resulting equilibrium 
coronal temperatures for the three observations are « 600 keV, 
« 2.2 keV, and « 1 .3 keV, respectively. When combined with 
their equilibrium optical depths, we find that the resulting 
Compton y-parameters range from 0.1-0.2. That is, Comp- 
tonization is only inducing a modest shift in the average pho- 
ton energy. However, given the fact that we do not have a good 
spectrum above 8 keV, there is not a strong enough lever arm 
on these Comptonization fits to say how free from systematic 
uncertainty these parameter values are. 

Likewise, detections of a (broad) Gaussian line peaking 
between 6.2-6.9 keV are marginal. The upper limits of the 
line strengths found here (in terms of integrated photon flux) 
are comparable to the estimated values from our prior RXTE 
observations, after taking into account the fact that RXTE 



also integrates line emission from the galactic ridge and 
likely had remainin g systematic Fe line residuals at the w 1 % 
level ( Nowak et al.l 12008 ). That is, these line strengths are 
marginally consistent with the RXTE results. Furthermore, 
the line fitted here may be partly substituting for general spec- 
tral hardening at > 6 ke V. It is clear, however, that any real line 
at the levels suggested by the RXTE observations cannot be 
narrow; any such line must be broad and have a peak residual 
relative to the continuum of < 2%. We illustrate this in Fig.|6] 
wherein we have refit the spectra excluding the Fe line, and 
combine all three observations in the plot — but not the fit — 
to show the Fe region residuals. 

The most salient feature of these fits is that a model with 
a single normalization constant is both successful, and yields 
a fairly small normalization of « 1.9 x 10""*. The eqpair 
model normalization is proportional to the disk flux that 
would be observed in the absence of Comptonization, and 
is the same as the normalization of the diskpn model 
jGierliriski et alJI 19991) . namely 



-^eqp ~ 



M 



D 
Ikpc 



/T COSi 



(1) 



where M is the central object mass, D is the distance, / is the 
disk inclination, and /c is the color-correction factor of the 
disk (see also Coppi 1999). For a mass of 3M0, a distance of 
lOkpc, and an inclination of / = 75°, we expect a normaliza- 
tion of w 3 X 10"^, a factor of 15 times higher than observed. 
(The discrepancy only increases for a more face on geometry; 
a more highly inclined geometry is unlikely given the lack 
of eclipses.) The fitted normalization can be recovered if we 
increase the distance to 40 kpc, increase the color-correction 
factor for the seed input photons to 3.3, or we decrease the 
characteristic emission radius by a factor of four (The emit- 
ting area dependence is responsible for the term in the 
normalization.) 

The diskpn model, which provides the seed photons 
for the eqpair model, has its peak emission radius at « 
IQGM/c^. If this radius were shrunk to 2.5GM/c^, the fitted 
normalization can be recovered with M = 3Mq, D = 10 kpc, 
and / = 75°. This is what we essentially now do by consider- 
ing model fits with the kerrbb model. The kerrbb model 
(iLi et al.ll2005h uses a relativistic model of the disk structure, 
and includes Doppler beaming and gravitational redshift in 
calculation of the emergent spectrum. The spin of the black 
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Fig. 7. — The joint fit of a simplCSkerrbb model to the three spectra 
assuming a black hole mass and distance of 16Mq and 22kpc. The spectra 
share a common interstellar absorption column, black hole mass and spin, 
disk inclination, and a fixed color correction factor of /c = 1-7. They have 
independent disk accretion rates and simpl parameters. Black/grey spectra 
are the combined XIS and 3 spectra, and blue/light blue are the XIS 1 
spectra. Lightly shaded data were not included in the fit. The top residuals 
panel shows the ratio residuals for obsID 405057010, the middle residuals 
panel shows the ratio residuals for obsID 405057020, and the bottom shows 
the ratio residuals for obsID 405057030. 

hole becomes a fit parameter. Increasing the spin from 
(i.e., a Schwarzschild black hole) has two major effects: it 
increases the accretion efficiency of the disk, and decreases 
the disk inner radius. These two effects in combination tend 
to increase the disk temperature while decreasing its normal- 
ization (i.e., emitting area or, equivalently, inner radius). Both 
of these effects are suggested in the 4U 1957+1 1 spectra. 

Here we adopt the exact same set of parameters as we 
considered fo r our prior studies of the RXTE spectra of 
4U 1957+11 (iNowak et al.ll2008h . We assume a disk incli- 
nation of i = 75°, the theoretically preferred color-correction 
value of /c = 1.7 (|Davis et al. 2005, 2006) for the disk spec- 
trum, and choose two sets of parameters: 3Mq and lOkpc, 
and 16Mq and 22kpc. We discuss these latter two choice 
more fully in Section |4] We further modify these fits by al- 
lowing for a spectral hardening which we model with simpl 
dSteiner et al. 2009), which is a convolution model which 
mimics some aspects of optically thin ComptonizatiorQ. The 
two parameters of this convolution model are the slope of the 
powerlaw into which disk photons are transferred, and the 
fraction of input photons transferred into the hard tail. 

As for the Comptonization models, we fit all three observa- 
tions simultaneously, include a broad Gaussian line in exactly 
the same manner, and tie some of the model parameters to- 
gether Specifically, we constrain each observation to be fit 
with the same neutral column, line energy and width (fixed 
to 0.3 keV), and black hole spin, but let the disk accretion 
rate, the simpl model parameters, and the line normaliza- 
tion be free. Results for these kerrbb model fits are pre- 

^ simpl, however, should not be confused with a Comptonization model. 
Aside from lacking any physical parameters, e.g., a coronal optical depth, a 
coronal temperature or compactness, and any reference to coronal geometry, 
it also cannot mimic spectral shapes associated with optically thick Comp- 
tonization. It should be regarded as purely a phenomenological model, albeit 
one that is more useful than a powerlaw since it does not extend its spectrum 
to low energies in the same manner as a steep powerlaw can. 



sented in Tableland Fig.|2] The resulting fits are good, with 
the 3M0, lOkpc model yielding xl/i' = 2063/1014 and the 
16Mq, 22kpc model yielding xV^ = 1830/1014. The latter 
model is very comparable to our best Comptonization fit. Ad- 
ditionally, this latter, more successful, model has line param- 
eters very similar to those from the Comptonization fits. We 
also find that only « 7- 15% of the disk photons are trans- 
ferred into a hard tail. The fitted neutral column is slightly 
smaller than for the Comptonization fits, but only by < 10%. 

Just as the Comptonization models could fit all three obser- 
vations with a constant normalization (equivalent to a constant 
inner disk radius), so too can the relativistic disk models fit all 
three observations with a constant spin parameter (also equiv- 
alent to a constant inner disk radius). The main variable pa- 
rameter among the observations here becomes the disk accre- 
tion rate. For both sets of mass and distance assumptions, the 
fitted spin parameter is high: nearly 0.9 for the 3Mq, lOkpc 
case, and nearly maximally spinning in the other We discuss 
the scaling relations among the disk parameters in the next 
section. 

4. SCALING RELATIONS 

As shown above, spectral fits to the data are somewhat de- 
generate. We have found a successful set of relativistic disk 
model fits assuming a mass of 3Mq and a distance of lOkpc, 
and also when assuming a mass of I6M0 and a distance of 
22kpc. The first set of parameters was chosen as being the 
lowest mass/closest distance black hole that is likely to be 
consistent with 4U 1957+11 persistently remaining in a soft 
state. The second set of paramete rs was chosen to be consis- 
tent with our prior RXTE studies (iNowak et alJl2008l) . which 
in turn were found via the scaling relations described below. 

As the spectra are dominated by the seed photons from a 
disk for both the disk and Comptonization models, the flux, 
F, scales as: 



4 /Tcoi 



(2) 



where D is the distance to the source, R is the inner radius of 
the disk, M is the black hole mass, Tgff is the disk's peak effec- 
tive temperature, Tcoi is the disk's peak color temperature, and 
/c is the color correction factor. The effective temperature to 
the fourth power scales as the fractional Eddington luminos- 
ity, L/LEdd, multiplied by M/R^ oc M"'. Thus, if we make 
the assumption that the lowest flux RXTE observations de- 
scribed by Nowak et al. ( 2008) are indeed at a fixed fraction 
of Eddington luminosity, e.g., 3%, then we will have: 



Tcol 
fc 



ocM~ 



(3) 



The color temperature is fixed by observation, and L/Lndd « 
M/M (where M is the accretion rate) is fixed by assumption, 
leading to the first two of our scaling relations: 



Ma/4 , M(x/^ 



(4) 



Just as the color temperature is fixed by observations, so too 
is the flux. Given that the disk inner radius scales as mass, and 
using the fact that the color temperature and flux are fixed by 
observation, we have that M^/D^ f* is fixed. Using the above 
scaling relations, we then obtain: 



(5) 
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We have previously used these scaling relations to find 
spectral fit s for a black hole m ass of 16Mq and a distance of 
« 22kpc (Nowak et al."2008'). Specifically, we first explored 
fits to RXTE spectra where we fixed the mass and distance to 
3Mq and lOkpc. We found that the best fits clustered around 
a spectral hardening factor of « 1 . 1 . We then fixed the spec- 
tral hardening factor to /c = 1.7 and the mass to 16Mq. We 
allowed the distance to be a free parameter, and found that its 
best fit values clustered around 22 kpc, i.e., roughly consistent 
with the scaling relations above. 

Here we perform a similar exercise for the Suzaku spectra. 
As in Section [3] we fit all three spectra simultaneously with 
the simpl(g)kerrbb model wherein we only allow the ac- 
cretion rate and simpl parameters to vary among the three 
spectra. We consider two sets of fits, one with mass of 3 M© 
and distance 10 kpc, and another with 16 Mq and distance 
22 kpc. For both we allow the black hole spin to be a free 
parameter, and we create grids of fits for fixed spectral hard- 
ening factors ranging over /c = 1 • 1-3.7. The plot of fitted spin, 
a*, vs. spectral hardening factor, /c, is presented in Fig. [8] 

As for our prior RXTE observations, we find that fits using 
3Mq, 10 kpc, and /c = 1.1 or 16Mq, 22 kpc, and /c = 1.7 
are nearly equally good and are of comparable to our best 
Comptonization fits. We also see that all fits using 16Mq, 
and 22 kpc require large spin parameters. The fits using 3Mq 
and 10 kpc only achieve low spin parameters for large color 
correction factors, /c > 3. (These conclusions are unaltered 
if uniform 1.35% systematic errors are included in the fits.) 

5. DISCUSSION 

There are four salient features of these Suzaku observations 
of 4U 1957+1 1 that echo the results that we obtained with our 
prior Chandra, XMM-Newton, and RXTE studies. The peak 
disk temperatures are high, the disk normalizations (or equiv- 
alently, emitting areas or inner disk radii) are small, all ob- 
servations can be fit with the same inner disk radius, and the 




lis 2 2.5 ' ' 3 3.5 



^0 I c' I Err 

Fig. 8. — Results of fitting a simpl{55kerrbb to the three Suzaku ob- 
servations of 4U 1957+11. The top panel shows the fitted spin vs. spectral 
hardening factor The bottom panel shows the reduced x" value of the fits. 
The blue (thicker) lines represent models with the mass and distance fixed to 
16Mq and 22 kpc, respectively, and orange lines represent models with the 
mass and distance fixed to 3 Mq and 10 kpc, respectively. 

spectra are remarkably simple. The only significant devia- 
tions from a simple continuum, aside from spectral hardening 
and a possible weak, broad Fe line at high energy, are « 0.7- 
0.9 keV features associated with absorption by the ISM, pre- 
dominantly due to the Fe L and Ne edges. 

We show a closeup of this region in Fig. |9] The spec- 
tral feature seen at w 0.55 keV, not described by our fitted 
models, can in fact be seen in the flux corrected spectra of 
any monochromatic, narrow line spectrum simulated using 
the Suzaku response functions. Small calibration errors in 
the amplitude of this feature potentially could lead to a sys- 
tematic underestimate or overestimate of the predicted counts 
near 0.55 keV. We thus do not believe that this feature is in- 
dicative of physical complexity in the underlying 4U 1957+1 1 
spectra. Likewise, it is unclear whether the w 5% residuals 
seen near the 0.7 keV Fe L edge represent complexity in the 



TABLE 3 

constant*TBnew*simpl(g) (kerrbb+gaussian) FITS TO 4U 1957+1 1 SPECTRA 



ObsID 



A'h 



M 



^xis3 



x7dof 





(10^' cm-2) (Jc/GM^) 


(10" gs-') 






(10^) 






405057010 
405057020 
405057030 


n 1 Qs+o 002 n sofis+o oou 
U.185_(,om 0.8965^0012 


0.180:«;»«i 
0.131!«;»«1 


1 1+-' 
1 ,+0.1 
^■^-0.0 

4.0+?;0 


000+0 001 

037+0.004 

015+0.006 
u.uij 0.008 


n 9+0'' 

0+<' l 

0+0.2 
"■"-0.0 


1 019+0 002 

1 nvR+o.oos 

1 067+0.005 
i.uo;_ooo4 


1.020!};-5][;| 2062.5/1014 
1 ns^+o'oo3 

1 n';7+o.oo4 

l.UJ/_ooo5 



405057010 
405057020 
405057030 



0.186+ 



n RQ<;7+0.0016 
0.895 /_o^0016 



n loi+o.ooi 

"•^^^-0.002 

isi+o.ooi 

"■^-^'-0.002 
I 1 (:+0.002 
"•^^"-0.00 1 



1.1+ 



1+2.9 
-0.0 

^■^-0.0 

4.0+00 



000+0 002 
"■"""-0.000 

0.032!«^005 



0.009 



-0.006 
1+0.008 
-0.005 



1.0+ 



1+0.5 
-0.5 

0.0+03 



1.021+ 



l+fl.003 
-0.004 
1 077+0.004 
^■""-0.004 
1 067+<'006 



1.021+ 



I +0.003 
^-0.004 

1 na 1+0.004 

^■^°i-0.004 

, 056+0.006 



1189.8/1014 



405057010 
405057020 
405057030 



17J+O.002 
"■" ^-0.002 



nqqq 1+0 0002 
"■''^^-0.0001 



n 7SS+0.007 
"■"^^'^-o.oos 
224+0.005 
"■^^^-0.004 
i7f:+0.003 
"■^ '"-0.004 



1 4+0.1 
^■^-0.1 

1 1+0- ' 
^■^-0.0 
o .+0.6 
-'■^-0.4 



072+0.009 

"■"'■^-0.007 
] 50+0 007 
u.l3U_ooo7 
1 27+0.032 
"■^^'-0.022 



2 3-^-^ 

^■^-0.5 
,+0.4 
^-0.1 
7+0.6 
"■'-0.5 



0.1+ 



1 021+0.003 
'■"^ -0.002 

1 nsn+ooo3 

1 069+0.005 
i.uoy 0.004 



1 021+*oo3 

i.o^i_0oo2 

, nss+0.003 

'■"'^^-0.003 

1.058^'!!04 



1830.0/1014 



405057010 
405057020 
405057030 



171+0.002 
"■" ^-0.002 



n 9Q9:!+0.0002 
u.777j_(j 0002 



n 7 8 7+0.009 
"■^° -0.007 
n 219+0.007 
■ -0.006 
n 170+0.004 
-^-0.004 



1 4+0.2 

1 1+0.' 

^■^-0.0 
q 7+0.7 
•'■ '-0.3 



0.068!fl:Oio 

"■142!":™ 

O.I28!0:026 



2 1+0.6 
^■^-0.5 

02+0.5 
"■^-0.2 
9+0.6 
"■^-0.6 



1 n77+O.004 
i.Ui^_Ooo3 

1 nRn+o.oo4 

l.USU 003 
1 nfiR+0.007 
^■Oo»-0.005 



1 077+0.003 
l.uzz_ooo4 

1 n»q+0.004 

'■'''^•'-0.004 

1 057+<''00« 

i.oj/_ooo6 



1084.3/1014 



Note. — Errors are 90% confidence level for one interesting parameter A'h is the neutral column; a* is the black hole dimensionless 
spin; M is the disk accretion rate; /c, the disk color-correction factor, has been fixed to 1.7; F and /scat are the simpl parameters 
(powerlaw slope, constrained to lie between 1.1^, and fraction of scattered continuum); is the Gaussian line normalization in units 
of photons cm"^ s"' ; and QisN are the normalization constants for the XIS and 3 detectors, relative to the XIS 1 detector. The top set 
of parameters are for a black hole mass of 3Mq and a distance of lOkpc, while the bottom set are for a mass and distance of I6M0 and 
22 kpc. Within each set of parameters, the first group ai'e for data that include only statistical errors, while the second group are for data 
that also include 1.35% systematic en'ors. 
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Fig. 9. — The same fit as seen in Fig. [5] but sliowing tlie XIS 0, 1, and 3 
spectra combined, and the model applied down to 0.5 keV. 

physical spectrum, as opposed to uncertainty in the response. 
Aside from the Fe L edges, Chandra and XMM-Newton spec- 
tra di d not show any spe ctral features near the 0.7 keV band- 
pass (iNowak et al.ll20()8l) . 

The only clear spectral features are the Ne edge/Ne IX 
absorption near 0.9 keV. Similarly simple disk spectra were 
implied by Suzaku observati ons of LM C X-3 in its soft, 
dis k dominated s tate (K ubo taetalJ 1201 0). As pointed out 
by iKubota et al.l (1201 Oh . such spectral simplicity is in fact 
surprising if modern disk atmosphere models apply to these 
sources. These disk atmosphere models predi ct spectral struc- 
ture with ±5% amplitude near metal edges (Davi s & Hubenvl 
11)06) that is seen neither here nor in Suzaku observations of 
LMC X-3. That the simplest, phenomenological disk mod- 
els lacking such features describe the spectra extremely well 
(i.e., to within < 5% 0.7-0.8 keV, and to < 2% above 0.8 keV) 
is quite remarkable, and may indicate that additional physics 
needs to be i ncorporated into the disk atmosphere models 
(IKubota et al. 2010). 

We also have confirmed the results of our prior studies in 
that the best fits assuming a small mass and distance {3Mq, 
lOkpc) require a black hole spin a* > 0.9. The required spin 
increases if we instead assume 16Mq and 22kpc. Altering 
the assumed disk inclination does not change this situation. 
Higher disk inclinations are ruled out by the lack of system 
eclipses. Lower disk inclinations allow for more observed 
gravitational redshifting, and less relativistic beaming, of the 
spectrum, and thus require larger spectral correction factors. 
(I.e., inclination effects more than just the projected area of 
the disk.) As an example, we took the 16 Mq and 22kpc 
kerrbb model, froze the inclination to 45° and let both the 
black hole spin and hardness correction factor become free 
parameters. A solution with comparable to the best fits 
discussed above was found, however, it required maximally 
negative spin of a* ss -1 and an extremely large color cor- 
rection factor of « 8. Thus, high spin/high inclination solu- 
tions are the most "straightforward" solutions that we have 
been able to obtain, and generally have had the lowest color- 
corr ection factors of th e models that we have explored (see 
also lNowak etallllOOSh . 

These results are in fact confirmed by the Comptoniza- 
tion model fits: the high disk temperatures (1. 31-1.48 keV) 
and especially the low normalization — 15 times lower than 
expected for a 3Mq Schwarzschild black hole at lOkpc — 
could be indicative of a rapidly spinning black hole. As dis- 
cussed above, the disk temperature in the diskpn seed pho- 
tons used in the Comptonization model peaks at w 10 GM/c^. 



If one were to reduce this radius to w 2.5 GM/c^, i.e., change 
the emission profile to more closely match that of a rapidly 
spinning hole, then the Comptonization model normalization 
would be consistent with a 3Mq black hole at only lOkpc. 

The Comptonization models offer two possibilities for fit- 
ting the 4U 1957H-1 1 spectrum with a low spin. The first, per- 
haps less likely, possibility is to increase the color-correction 
factor via a low temperature, optically thick corona. We have 
already seen that the lowest flux observations are fit with such 
coronae, albeit with a Compton y-parameter of ss 0. 1 . Instead, 
we are envisioning a scenario wherein y w /c « 3 - the color 
correction value required to fit a 3Mq, lOkpc, Schwarzschild 
black hole. This would require an optical depth of 

where kT^ is the coronal temperature. This is a coronal regime 
typically not explored in spectral models. The concept would 
be that the low flux, equilibrium solution of the disk was such 
a cool, optically thick atmosphere which at higher flux (i.e., 
the 15% of the time requiring a hard tail) expanded into a hot, 
optically thin corona. 

Aside from lacking a theoretical description of such a sce- 
nario, the fitted constant emitting area at low flux implies that 
the 4U 1957H-1 1 system must always return to the same equi- 
librium coronal solution at low flux. It is not clear why this 
should be the case. The second possibility offered by the 
Comptonization solutions is simply to increase the distance to 
4U 1957H-1 1, while keeping the mass small. The scaling rela- 
tions discussed in Section|4]only apply under the assumption 
that the lowest flux observations are a fixed fraction of the Ed- 
dington luminosity (e.g., 2%LEdd)- If we relax this assump- 
tion, then eq.[T]suggests that if the 4U 1957H-1 1 distance were 
40kpc, then it could be a 3Mq Schwarzschild black hole. 
This would imply, however, that the faintest 4U 1957H-11 
spectra are > 30%LEdd, while the brightest observed RXTE 
spectra are > Ledd- This would be rather unusual, espe- 
cially given the lack of strong variability found in the RXTE 
spectra: the root mean square variability is typically < 3%, 
and at its strongest was < 13% in the 2-60 keV bandpass 
over w 10-2-2 x 10^ Hz ( iNowak & Wilm s 1999; Nowak etiD 
l2008h . 

It therefore seems likely that 4U 1957H-11 is a system that 
has a combination of both high spin and large distance, which 
puts this X-ray binary in the Galaxy's halo. It is interesting 
to speculate whether these two characteristics are related to 
one another. First, we note the obvious fact that in order for 
4U 1957H-11 to be a X-ray binary it must still be accreting 
and therefore be a comparatively young system. Given the 
low star formation rate in the Galaxy's halo, it is rather un- 
likely, therefore, that 4U 1957H-11 was formed in the halo. 
Instead 4U 1957H-1 1 likely originates in the Galactic disk and 
then migrated into the halo. Several possible formation sce- 
narios have been discussed in the literature, most recently as- 
sociated with the discovery of young hypervelocity stars in 
the Galaxy's halo. For example. iPeretsI f2009, and references 
therein) discuss a scenario where hypervelocity binary sys- 
tems are created through the interaction of a hierarchical triple 

* Note that we are ignoring issues of the unknown inclination. Choosing 
a more face-on inclination only exacerbates the normalization problem (i.e., 
we would have expected an even larger normalization), and would introduce 
more gravitational redshifting to the observed temperature, increasing the re- 
quirement for rapid spin in the spectral fits. 
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star system with the supermassive black hole in the center 
of the Galaxy. Since hierarchical triples are very common 
among multiple star systems, such a scenario is likely. With a 
typical ejection velocity of ISOOkms"', it would take this sys- 
tem a few 10^ years to reach its current position, well within 
the lifeti me of typical evolutio n scenarios for X-ray binary 
systems (Ivan den Heuve]|[T976l) . (At such a velocity, and as- 
suming a 22kpc distance, the proper motion of 4U 1957+11 
likely would have amount ed to < 0.3'^ since t he discovery of 
the optical counterpart bv Margo n et al.|[T97 8.) Such a sce- 
nario could also explain a large spin for the system, if during 
the ejection angular momentum is transferred onto the black 
hole. 

A high black hole spin for 4U 1 957+ 1 1 is a consistent inter- 



pretation of the results presented here; however, such a result 
is dependent upon the unknown mass, distance, and inclina- 
tion of the system. Given the nature of 4U 1 957+ 1 1 as perhaps 
the simplest, cleanest example of a BHC soft state, and per- 
haps the most rapidly spinning black hole, observations to in- 
dependently determine this system's parameters are urgently 
needed. 
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